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Abstract —Due to the increasing amount of electricity gen¬ 
erated from renewable sources, uncertainty in power system 
operation will grow. This has implications for tools such as 
Optimal Power Flow (OPF), an optimization problem widely 
used in power system operations and planning, which should 
be adjusted to account for this uncertainty. One way to handle 
the uncertainty is to formulate a Chance Constrained OPF 
(CC-OPF) which limits the probability of constraint violation 
to a predefined value. However, existing CC-OPF formulations 
and solutions are not immune to drawbacks. On one hand, 
they only consider affine policies for generation control, which 
are not always realistic and may be sub-optimal. On the other 
hand, the standard CC-OPF formulations do not distinguish 
between large and small violations, although those might carry 
significantly different risk. In this paper, we introduce the 
Weighted CC-OPF (WCC-OPF) that can handle general control 
policies while preserving convexity and allowing for efficient 
computation. The weighted chance constraints account for the 
size of violations through a weighting function, which assigns a 
higher risk to a higher overloads. We prove that the problem 
remains convex for any convex weighting function, and for very 
general generation control policies. In a case study, we compare 
the performance of the new WCC-OPF and the standard CC- 
OPF and demonstrate that WCC-OPF effectively reduces the 
number of severe overloads. Furthermore, we compare an 
affine generation control policy with a more general policy, 
and show that the additional fiexibility allow for a lower cost 
while maintaining the same level of risk. 

1. Introduction 

Optimal Power Flow (OPF) is an optimization problem 
widely used in power system operational and expansion 
planning. In this paper we will address the application of 
OPF in operational planning. The objective is to minimize 
cost of operation, while avoiding violations of technical 
constraints such as limits on transmission or generation 
capacity. However, with increasing penetration of electricity 
generation from renewable sources, power systems operators 
face higher operational uncertainty. To avoid violations of 
technical constraints, it is increasingly important to account 
for these uncertainties during the operational planning phase, 
i.e., within the OPF problem. 

One way to handle uncertainty is to limit the probability 
of constraint violation by formulating a chance constrained 
OPF (CC-OPF). In [1], a joint chance constrained problem 
was formulated, and solved using the scenario approach, a 

L. Roald and G. Andersson are with the Power Systems Labora¬ 
tory, Department of Electrical Engineering, ETH Zurich, Switzerland 
{roaldIandersson}@eeh.ee.ethz.ch 

S. Misra and M. Chertkov, are with the Center for Nonlinear Studies and 
Theoretical Division T-4 of Los Alamos National Laboratory, Los Alamos, 
NM 87544. M. Chertkov is also affiliated with the New Mexico Consortium, 
Los Alamos, NM 87544. {sidhant | chertkov}@lanl. gov 


solution method based on samples. Both [2], [3] formulated 
the CC-OPF with separate chance constraints, and obtained 
an analytical reformulation assuming normally distributed 
wind power fluctuations. Convexity of the CC-OPF with 
respect to the characteristics of controllable generation was 
a key observation of [2] which also allowed for guarantees 
of exactness and lead to efficient implementation for large 
networks (thousands of nodes). 

While previous versions of CC-OPF have shown the ability 
to limit the probability of violations, they have some draw¬ 
backs. First, current CC-OPF formulations only consider 
affine policies for generation control, which is not necessarily 
realistic for large wind power deviations and may lead to 
sub-optimal results. Second, standard chance constraints do 
not distinguish between large and small violations, although 
those might carry signiflcantly different risk. This paper ad¬ 
dresses both drawbacks, thus extending existing formulations 
in two ways. 

First, we suggest more general policies for generator 
control. An affine policy is a good approximation of the 
existing implementations of automatic generation control 
(AGC) for small deviations. However, a more general policy 
allows us to model additional control actions, such as manual 
activation of reserves during large wind power deviations. It 
is thus possible to obtain a cheaper, more realistic solution 
while maintaining the same level of risk. To the best of our 
knowledge, no previous chance constrained OPF aimed at 
reducing the cost in a similar way. 

Second, the proposed extension accounts for both the 
probability and the size of constraint violations. A standard 
chance constraint only considers the probability of violations, 
and does not account for the risk related to the actual 
constraint violations. To remedy this, we deflne a weighting 
function for the overload, and calculate the risk as the product 
between the weighting function and the probability distribu¬ 
tion of the overload. This Weighted CC-OPF (WCC-OPF) 
is convex with respect to parameters of the generation, thus 
allowing exactness and tractability in the spirit of [2]. While 
the set-up allows for general convex weighting functions, 
the considered examples are chosen to be applicable to the 
power system problem. The examples are closely related to 
the severity functions applied in risk-based OPF (e.g., [4], 
[5], [6]), where they are used to account for the severity 
of post-contingency overloads. Such a severity function was 
combined with a CC-OPF in [6]. However, [6] limited the 
probability of high risk situations, and did not consider 
the average risk due to constraint violation. The weighted 
chance constraint suggested here is somewhat similar to a 



constraint on the conditional-value-at-risk (CVaR), a risk 
measure borrowed from finance, which was introduced by 
[7] and applied to the OPF problem in [8]. However, the 
weighted chance constraint suggested in this paper is more 
general and adapted to better fit the power system problem. 

The remainder of the paper is organized as follows. 
Section II introduces the original CC-OPF. In Section III, the 
weighted chance constraint is introduced, along with a proof 
of convexity. Section IV provides some relevant examples for 
practical weighting functions with affine and piecewise affine 
generation control. Section V illustrates the performance of 
the method based on the IEEE RTS96 test system. Section 
VI summarizes and concludes the paper. 

A. Notations 


only to the total deviation in wind generation output, as in 
[2], all columns of a are identical, i.e., cx(.^i) = 

3) Power Flows: The power fiows pij on each line are 
computed according to the standard so-called DC approxi¬ 
mation [9] 


Pij =M(ij,.)(p-aw + v + w-d),Vijg£. (2) 

The matrix M G relates the line fiows to the 

nodal power injections, which are expressed as the sum of 
generation p — cxu, wind power production vFuj and demand 
—d. M is defined as 


M = Bf 


(Btus) ^ 

0 


0 

0 


(3) 


We denote vectors by lower case letters p, cc. The compo¬ 
nents of the vectors are denoted by using subscripts, i.e, the 
ith component of p is denoted by pi . Matrices are denoted by 
upper/lower bold case letters, a, M, and denote 

the row and column of a, respectively. Index i refers to 
generators, while index ij refers to lines. 

II. Chance constrained Optimal Power Flow 

In this section, we review the existing literature regarding 
the modeling of the OPF problem with fluctuating in-feeds 
as a CC-OPF. Our description is based on the formulation 
in [2]. Let (G = (V, f) represent the graph of the power 
network, where V is the set of nodes with \V\= m and £ is 
the set of edges/lines of the system with \£\ = n. The set of 
wind generators is denoted by W C V. Wind is considered 
as the only source of fluctuations in this paper, although the 
formulation can easily be extended to handle fiuctuations 
from other sources such as solar PV or load. The set of 
non-wind generators is denoted by ^ C V, and are assumed 
to be controllable within their limits. To simplify notation, 
we assume that there is one controllable generator p, one 
wind generator w and one demand d per node, such that 
\Q\ = |W| = |V| = m. Nodes without generation or load 
can be handled by setting the respective entries to zero. 

1) Wind power in-feeds: The wind in-feeds w are mod¬ 

eled as the sum of the forecasted electricity production from 
wind, given by v = and a zero mean fiuctuating 

component uj = w — v. 

2) Generation control: Since secure operation of the 
power system requires balance between produced and con¬ 
sumed power at all times, any deviation uj in the wind power 
production must be balanced by an adjustment in the control¬ 
lable generation. In the CC-OPF [2], these adjustments are 
modelled through an affine policy, refiecting the automatic 
generation control which is establishing balance within tens 
of second to a few minutes [9]: 

p{uj) = p-aoj , (1) 

where p is the vector of scheduled generation, and a G 
j^mxm ^ matrix with elements describing the response of a 
generator to a wind fiuctuation. The elements of each column 
in a sum to one.If we assume that the generators respond 


where Bf is the line susceptance matrix and Bi)us the bus 
susceptance matrix (without the last column and row) [1]. 
.) is the row of M related to the line (ij) G £. 

4) Optimization Problem: With the above modelling con¬ 
siderations, the CC-OPF is stated as follows: 

min CiPi (4) 

s.t. (5) 

y2pi-di+vi = o (6) 

iev 

= l Vjgw, p>0, a>0 (7) 

ieg 

Pu [Pi - < e* '^ieg (8) 

[pi < p™*"] < ei Vigg (9) 

Pu, -aw-d + t>+w)> Pij] < Cij 

( 10 ) 

Pw [M(,j,.)(p -auj-d + v + uj)< -pij] < €ij '^{ijjes 

( 11 ) 

The objective Q is to minimize generation cost, where c 
contains costs, i.e. bids, from the generators. Eq. <0 refiects 
the power balance during nominal (pre-planned) operation, 
while 0 accounts for deviations from nominal, where p 
and a are assumed positive. (The latter constraints is not 
theoretically significant, but refiects engineering practice). 
Since the generation output and the power fiows depend on 
the fiuctuations cc, the generation constraints 0, 0 and 
transmission constraints ([^, ([^ are formulated as chance 
constraints with accepted violation probability e. To obtain a 
tractable optimization problem, the chance constraints must 
be reformulated. Assuming that the wind fiuctuations follow 
a multivariate normal distribution, 0-0 can be restated 
as second-order cone constraints [2]. 

The CC-OPF formulation above accounts for uncertainty, 
but have some shortcomings. First, it assumes that the gen¬ 
eration control follows an affine policy. Second, the chance 
constraints treat all violations equally, irrespective of the size 
of the violation. In the following, we discuss a new way 
of formulating the chance constraints to address both those 
shortcomings. 




III. Weighted Chance Constraints 


To formulate a chance constraint that can handle general 
policies for generation control and at the same time differ¬ 
entiates between large and small violations, we propose a 
generalized weighted chance constraint of the form 

/ oo 

< e • (12) 

-OO 

Here, P{uj) is the multivariate distribution function of the 
fluctuations. The quantity y{uj) denotes the magnitude of 
overload and depends on whether we are considering a 
violation of the upper or lower limit on generation or power 
flow. In case of violation of the upper generation limit for 
generator i, we have 

yiuj)=piiw)-pr-, (13) 

and in the case of lower limit violation, we have 

y{u!) -Pi{u). (14) 

Similarly, for power flows we have that for line {i,j) G £, 
y{Lo) = Pij{ui) - (15) 

and 

y{‘^)=Pij''^ ( 16 ) 

where 

Pij{to) = (]V[(ij_.)(p(w) -d + v+u)). (17) 


Overload y 
Standard f(x)=x(y>0) 



Overload y 
Linear f(x)=y5^(y>0) 


Overload y 


Fig. 1. Examples of weighting functions for chance constraints: Standard 
(left), linear (middle) and quadratic (right). 


Step function). While unit step function assigns equal risk 
to all magnitudes of overloads, the risk increases linearly 
or quadratically with increasing overload for the linear or 
quadratic weight functions. 

During large wind power deviations, the power system 
operator might use additional controls, such as slower re¬ 
serves (e.g., spinning or tertiary reserves), to balance the 
system and control power flows. Since these controls have a 
signiflcant influence on the ability to handle disturbances, 
they should be included in the OFF to avoid unrealistic 
or sub-optimal solutions. A main advantage of (H) with a 
convex weighting function f{y) is therefore in the principal 
ability to handle more general control policies while still 
maintaining convexity, as shown below. 

A. Convexity of weighted chance constraints 

We state our result regarding the convexity of the 
weighted chance constraint deflned in ( p^ in Theorem 


Whenever we have ^ > 0, it indicates a violation of the 
limit, while ^ < 0 implies that we are in a safe operating 
region. The weighting function f{y{uj)), which is nonzero 
only if y > 0, describes the risk related to the overload. Note 
that the constraint ( p^ depends implicitly on the generation 
control policy p{uj). 

The weighting function f{y) can be implemented in many 
ways. For example, ( p^ is equivalent to a standard chance 
constraint if f{y) is the unit step function, i.e., 0 for ^ < 0 
and 1 for ^ > 0. However, the step function is not convex, 
which means that the standard chance constraint will not 
always be a convex constraint. However, as we will see in 
Theorem [T] the constraint is a convex constraint for 
general generation control policies and general probability 
distributions of the wind fluctuations whenever the weight 
function /(.) is convex. 

While the desire to maintain convexity is a motivation 
on its own, a convex weight function, which assigns a 
higher risk to a higher overload, is also natural from an 
engineering point of view. For example, in the context 
of so-called risk-based OFF, linear [4] and quadratic [5] 
weighting functions have been used to model the risk of post¬ 
contingency overloads. Assuming that the risk of overload 
would be similar in the case of overloads induced by wind 
power fluctuations, weighting functions similar to these used 
in [4], [5] could be considered. For comparison, these weight 
functions are plotted in Fig. along with the weight function 
corresponding to a standard chance constraint (which is a unit 


Theorem 1: Let p{uj) be any general generator control 
policy. Consider a general weighted chance constraint de¬ 
flned in ( p^ given by 

/ oo 

f{y{(^))P{(^)duJ < e (18) 

-OO 

Assume that the weighting function /(.) is a convex function 
of its argument. Then ( p^ is a convex constraint in the 
control policy p{uj). 

Proof: Referring to Eqs. (p3])-(p^ and Eq. ( p^ , we 
see that the overload y{uj) is always a linear function of 
the generator control policy p. Combined with the convexity 
of the weight function /(.), this implies that the function 
f{y{(jj)) is a convex function of p for every uj. The proof 
then follows from the fact that the expectation of a convex 
function is also convex. ■ 

Theorem shows that the weighted chance constraints 
are convex under general control policies, as long as the 
weighting function /(.) is convex. In addition, since the the 
network must obey total power balance at all times, we must 
also have the additional constraint on the generation policy 
that enforces this, given by 

yy^Pi{uj)+Vi+LOi-di = 0, Wlo. (19) 

i 

Note that for every oo, the constraint in ( p^ is simply a 
linear constraint. This means that the control policy p can 













be a variable within a convex WCC-OPF problem. However 
for practical implementation of the weighted chance con¬ 
straints within a standard convex optimization program, it is 
necessary to represent the control policy with finitely many 
parameters. It is straight forward to see that the convexity of 
© and the linearity of is preserved as long as we use 
any linear-in-parameter representation of the control policy 
p given by 

K 

p{uj) = p-'^akgk{oj), (20) 

*=1 

along with the constraints 

p^v — d = 0, ( 21 ) 

K 

+ w = 0, (22) 

fc=i 

where ak G and gk{(^) ^ The original affine policy 
O can be obtained as a special case of ( [2Q| by setting K = 
m and gk{uj) = cuk. 

In the above, ^/c(cc) can be any general function of uj. 
The expression for the control policy in ( [20| ) has extensive 
representation potential. For example, we can represent most 
general control policies p{uj) (e.g., continuous or piece-wise 
continuous) approximately in the form of ( [20| ) by using a 
piece-wise constant representation of p{uj), i.e., 

K 

p{u)) =p-'^akxi‘^ ^ Sk), (23) 

k=l 

where Sk are sets that form a partition of the domain of 
UJ and ak is the value of p on Sk- Increasing the number 
of terms K in the representation allows us to approximate 
general control policies with higher fidelity, but the number 
of variables to optimize over in the WCC-OPF also increases. 
Thus, a trade-off between fidelity and computational effort 
must be made. 

IV. Weighted chance constraints with example 

WEIGHT FUNCTIONS AND CONTROL POLICIES FOR 
NORMALLY DISTRIBUTED FLUCTUATIONS 

In this section, we derive the expressions for the weighted 
chance constraints for linear and quadratic weight functions 
and for affine and piece-wise affine control policies. The 
expressions assume that the wind fiuctuation uj is distributed 
as a multivariate normal distribution. For completeness, we 
also contrast these expressions with the original chance 
constraint in our current notation. 


be handled similarly. Substituting p{uj) = p — oluj in ( p^j ), 
we get that the generator overloads y'f are distributed as 
Gaussian random variables with the average overloads 
and variances (crj^)^ given by 


Pi = Pi-Pi 


(O' = • (25) 


Similarly, using (15), we get that the line overload yfj is 
Gaussian with 


l^ij = M(ijv) {p-d + p.)- Pij , (26) 

(<t“ )' = . (27) 

Similar expressions can be obtained for the lower limits 
yl, ylj using (14) and (16). We observe that y is distributed 


as a Gaussian random variable in all cases, since they are 
linear combinations of the random Gaussian uj. 

Since ^ is a Gaussian random variable, the LHS of ( [24| ) 
is simply the expectation of a truncated Gaussian random 
variable, and can be expressed as 




1-^ 


ZjS 

(j 


+ 




D 2 \ a ) 


< e , (28) 


where ^{x) is the cumulative distribution function of the 
standard Gaussian. The weighted chance constraints for 
generator and line overloads are obtained by using the cor¬ 


responding means and variances from Eqs. ([25])-([27]) above 
in ([28l). 


B. Quadratic weight function and affine policy 
Using f{y) = y'^xiv > 0) in we get 

/ OO POO 

y'^xiv > 0)P{uj)du! = y^ P{y)dy < e. (29) 

-OO J 0 


The LHS of ( [29] ) is the second moment of a truncated normal 
and can be rewritten as 


As for the case of a linear weight function, the weighted 
chance constraints for the quadratic weight function can 
be obtained by substituting the means and variances from 
Eqs. in 


C. Standard chance constraints and affine policy 

Here we include the expressions for the standard chance 
constraints with f{y) = x{y> 0) for the sake of complete¬ 
ness. 


A. Linear weight function and affine policy 

For the linear weight function, substituting f{y) = 
yx{y > 0) in we get 


/ OO POO 

yx{y > 0)P{uj)du! = y P{y)dy < e. (24) 

-OO J 0 



x{y > o)P(w)c(w 


f 


P{y)dy < e 


which can be reformulated as 


/i + ^(1 ~ ^ 0: 


(31) 


(32) 


We derive the expressions for constraint violations of the 
upper generator and line limits. The lower limit case can 


and then the final form can be obtained by using Eqs. ( j^ - 
(27) as usual. 











D. Linear weight function and piece-wise affine policy 

In this section, we consider a piece-wise affine policy 
which can be expressed as 

{ P — OLUJ^ VL~ < 

p — OLUJ + 13^, ft > ft~^ > 0 (33) 

p — cxuj 3- Pf ^ ft < ft~ < 0 


Here, ft = "f2iev wind power deviation, and 

p~ represent additionally deployed reserves in case of 
large wind fluctuations, when ft is larger (or smaller) than a 
given threshold ft~^ (or ft~). To ensure a balanced system, 
we additionally enforce = ^i^g P~ = 0- 

Under a piece-wise affine policy, the linearly weighted 
chance constraint can be written as 


+ 

+ 


pOO pOO pOO 

/ yP{y)dy= / / yP{y\^)P{n)dyd^ 

Jo J — oo J0 

/ O pOO 

/ yP{y\n)P{n)dydn 

-oo J 0 

POO 

/ / yP{y\n)P{Sl)dydn 

Jq- Jo 

pOO pOO 

/ / yP{y\n)P{n)dydn. 

Jci+ Jo 


(34) 


(35) 


The random variable is normally distributed and its mean 
and variance depends on which of the three regions in ( [3^ 
the total wind deviation ft belongs to. For generator power, 
the conditional mean and variance are given by 

Air(o) -pr"+A(^) - n, 

(36) 


where 


pm = 


A”, 

0 , 

[Pt^ 


Q < 

ft- <ft< U+, 
ft > 0+. 


(37) 


E. Quadratic weight function and piece-wise affine policy 

Following the same derivation as in Section [rV-D[ we can 
obtain an expression analogous to ( [35] ) by simply replacing 
^ by . The rest of the derivation is identical with the only 
difference that we substitute the expression for the second 
moment of a truncated Gaussian from ( [^ . 

V. Case Study 

The WCC-OPF is illustrated in a case study based on the 
IEEE RTS96 system, with modifications to the base case 
[10] similar to these described in [3]. The line capacities 
are reduced to 80 % of the steady state capacity listed in 
[10], since only active power flows are considered. The 
generation costs used in the OFF are linear costs adopted 
from [11], and we also assume that the generation capacity is 
twice larger than those from the original RTS96 model [10]. 
The minimal generator output is set to 0 for all generators. 
Aggregations of wind power plants are located at bus 8 
and bus 15, with a forecasted in-feed of 125 and 175 MW, 
respectively. The wind power fluctuations are assumed to 
follow a multivariate normal distribution with zero mean and 
standard deviations of 9.4 and 13.1 MW, respectively, with 
a correlation coefficient p = 0.2. Finally, we assume that the 
generation policy only reacts to the overall wind fiuctuation 
ft, such that ([T]) reduces to 


p{uj) = p — aft . (41) 

Here, a G is a vector with one contribution factor for 
each generator. 

In the following, we investigate two cases. First, we 
assume an affine policy for generation control, and assess 
how the weighted chance constraints perform compared with 
standard chance constraints in terms of number and size of 
the constraint violations. Second, we investigate how a the 
more flexible piece-wise affine policy for generation control 
( [33] ) performs compared with the standard affine policy 0 
in the case of the linear WCC-OPF. 


and 

(a(,„)Sl)' . (38) 

Similarly, we can calculate the conditional means and 
variances of the line overloads as 

= M(i,..) \p-d + v + - pm 

+ (Ysi)”' .){/ - a)I]l) n (39) 

= M(y,,)(/ - a)E(/ - afMm 
- (YEl)“'(M(i,..){J-a)Sl)% (40) 

with P{ft) defined as in ( [JT] ). 

The final version of the linear weighted chance constraints 
is obtained by substituting the expression for the expectation 
of a truncated Gaussian from ^2^ in ( [35] ), and insert the 
means and variances from @ - Q. 


A. Comparison of CC-OPF and WCC-OPF 

We compare the standard CC-OPF with the linear and 
quadratic WCC-OPF assuming an affine policy for gen¬ 
eration control. We choose eij = 0.1, for the 

transmission line constraints in all settings discussed. For the 
generation constraints, we choose different for the different 
problems, with = 0.001, for the CC-OPF and the 

linear WCC-OPF and = 10“^, \/{i}^g for the quadratic 
WCC-OPF. 

The cost of the OFF solutions with the aforementioned 
Ci, eij is very similar for the three formulations, as shown in 
Table[I] The cost difference is less than 0.1%, with the chance 
constraint being slightly cheaper than the other two, and the 
quadratic weighting function being the most expensive. Since 
the cost of the different formulations is so similar, it is useful 
to compare their performance in terms of number and size 
of violations. 





TABLE I 

Cost for the CC-OPF and the WCC-OPFs 



CC-OPF 

Linear WCC-OPF 

Quadratic WCC-OPF 

Total cost [$] 

16546 

16547 (+0.01%) 

16562 (+0.1%) 


0.1 


Q. ro 

E o 2 

LU > Q. 



Transmission line 12 


L 

< 


L. 

b 


Standard 


Linear 


Quadratic 


To test the OFF solutions, we generate 10 000 samples 
from the multivariate distribution of the wind power fluc¬ 
tuations and compute the number and size of the resulting 
violations. For the generation constraints, the choice of a very 
small ei ensures that there are very few significant violations 
of the generation limits. However, since Ci > 0, the linear and 
quadratic WCC assign a small, but non-zero a to some gen¬ 
erators procuding at the limit. With symmetrically distributed 
wind power deviations with zero mean, the corresponding 
generators thus experience a violation probability of 50%, 
but with very small violations <0.1 MW. 

For the transmission line constraints, where e is larger, 
there are more violations. We only consider the active 
transmission lines constraints, since this is where the effect of 
the constraint formulation is visible. The active transmission 
constraints are also the only ones who experience a signifi¬ 
cant number of violations. In Fig.[^ the empirical probability 
Eg of violations larger than 0, 1, 2, 5 and 10 MW is plotted 
for three transmission lines, line 12 (top), line 23 (middle) 
and line 28 (bottom). For the standard CC, the probability 
of any violation > OMW is close to 0.1 for all three lines. 
The linear WCC lead to a higher violation probability for 
line 12 and 23, but lower violation probability for line 28. 
The quadratic WCC is more restrictive, with all violation 
probabilities below 0.1. Considering the size of the different 
violations, we see that line 23 has some violations > bMW. 
For this line, it is clear how the standard CC allows for more 
and larger violations than the linear and quadratic WCC. 
Particularly the quadratic WCC reduces the number of large 
overloads to almost zero. 

Based on the results, we see that the WCC are less 
restrictive than the standard CC for small violations, and Cg 
exceeds the probability of the corresponding CC. However, 
both the linear and quadratic WCC are more effective in 
reducing the probability of large overloads, as seen for line 
23. This demonstrates the purpose of enforcing a weighted 
chance constraint as opposed to a standard chance constraint. 

B. Comparison of CC-OPF and WCC-OPF 

We now compare the linear WCC-OPF with an affine 
policy 0 and a piecewise affine policy ( [^ for generation 
control. We use the same set-up as above, but choose = 
0.01, and = 0.01, for both policies. 

For the piecewise affine policy, we define the threshold for 
deploying additional reserves /3+, (3~ as ^2+ = = 

70 MW. 

The cost of the WCC-OPF solutions are shown in Table 
1^ The cost of the CC-OPF with piecewise affine generation 
control is lower than the policy with affine control. The 
reason for this cost decrease is due to more fiexible use of 
generation resources. Fig. shows the generation output of 




Fig. 2. Empirical violation probability Ce for transmission line 12 (upper 
part), line 23 (middle) and line 28 (bottom). The results for all three 
formulations (standard, linear and quadratic) are shown in each plot, and 
the color of the bar indicate the empirical probability of exceeding certain 
violation thresholds (0, 1, 2, 5 and 10 MW). 

TABLE II 

Cost for the affine and piecewise affine policies 



Affine 

Piecewise Affine 

Total cost [$] 

16546 

16569 (-0.2%) 



Wind fluctuation Q [MW] 


Fig. 3. Generation output of generator 7 as a function of the total wind 
fluctuation $2 for the affine and the piecewise affine policy. 

one generator as a function of the total wind fiuctuation ft 
for the affine (straight line) and the piecewise affine policy 
(straight line with jumps). We observe that both policies 
lead to similar set-points for the scheduled generation p, 
and that they hit the maximum and minimum generation 
bounds at the same time. However, the piecewise affine 
policy allows a steeper slope than the affine policy, due 
to the use of additional reserves /3~^ ^ /3~ when the total 
wind fiuctuation exceeds the thresholds 0+, . For the 

congested transmission lines 12, 23 and 28, a similar effect 
is observed. For large fluctuations, the additional reserves 
shifts the power fiows to lower values, thus lowering 

risk. 

VI. Conclusions 

In this paper, we introduced the WCC-OPF as an extension 
to previous CC-OPF formulations. We defined a new type of 
chance constraints, the weighted chance constraints, that can 
















































handle general control policies while preserving convexity 
and allowing for efficient computation. The weighted chance 
constraints account for the size of constraint violations 
through a weighting function, which assigns a higher risk 
to larger violations. 

The performance of the new WCC-OPF was demonstrated 
in a case study for the IEEE RTS96 system. It was shown that 
the WCC-OPF allows for a larger number of small violations, 
but reduces the number of severe overloads compared to 
the CC-OPF. Further, it was shown that a more flexible, 
piecewise affine policy for generation control allows for 
lower cost compared with the original affine policy, while 
maintaining the same risk level. 

In our computations, both transmission and generation 
constraints were assigned the same type of weighting func¬ 
tion. It would however be possible to use different weighting 
functions for generation and transmission, since the conse¬ 
quences of a violation is different for the two. This type 
of analysis, along with an extension of the standard chance 
constraints to piecewise affine policies, is part of future work. 

Convexity of the WCC-OPF formulation discovered in 
this manuscript suggests a path forward for scaling up the 
approach to larger systems through deployment and develop¬ 
ment, in the spirit of 12], of a variety of computational gad¬ 
gets of the optimization theory such as cutting plane, piece- 
wise linear approximations of convex constraints and others. 
Notice also that our method allows generalizations to account 
for non-Gaussian effects in fluctuations of renewables, while 
preserving convexity of the formulation. Future work is 
needed to explore possible synergy of this observation with 
complementary work on improving statistical modeling of re¬ 
newables 112] and further exploration of advanced sampling 
techniques 113], 114]. The possibility to formulate general 
control policies can also be exploited to flexibly model other 
power system controls, such as FACTS or HVDC, extending 
the current affine models applied in 115]. Finally, we plan 
combining collection of ideas and techniques proposed in 
the manuscript with recent ideas based on energy function 
116], 117] and monotone operator 118] approaches to attack 
more challenging CC- and WCC versions of the nonlinear 
AC-OPF formulations, thus accounting for risks beyond the 
DC-approximation explored so far. 
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